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ABSTRACT 


The radiation output from the NPS Pulserad 112A X-ray Generator has 
been calculated using the ITS code CYLTRAN. The results of this numerical 
simulation have been compared to the experimental measurements taken previously. 
These calculations show that the overall photon energy spectrum is independent of 
the electron source beam parameters such as beam radius and angular distribution. 
A previously unexplained measured dip in the radiation dose at the beam center line 
can be explained with an angular divergence of the electron beam from the cathode. 
The simulation is successful in explaining the pattern of the dose distribution, but 
the calculation is in general much smaller than the measured values. This inability 
to reproduce the magnitude of the dose pattern points out the need to measure the 
radiation dose each time the Pulserad 112A is used. Some suggestions for explaining 


this discrepancy are included. 
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I. INTRODUCTION 


The Model 112A Pulserad Pulsed X—Ray Generator installed at the Naval 
Postgraduate School was built by the Physics International Company. Installation 
in the new Flash X-ray Facility was completed in August of 1988. This system is a 
high power x—ray source. The radiation is produced in a three step process. The 
twelve stage Marx Generator is charged via an external power supply to the desired 
voltage. The Blumlein Pulse Forming line is then resonance charged from the Marx 
Generator. When fully charged the Blumlein discharges rapidly into the electron 
accelerator tube creating a large potential difference across the diode gap. The 
cathode consists of a stainless steel rod which readily emits electrons when the high 
voltage is applied. The electrons are accelerated across the anode—cathode spacing 
until they impact a tantalum target. The resulting bremsstrahlung process produces 


radiation. 


А. PREVIOUS EXPERIMENT ON RADIATION OUTPUT 

The use of the Model 112A Pulserad Pulsed X—Ray Generator in radiation 
effects study requires some knowledge of the characteristics of the output radiation. 
The experiment to measure the magnitude of the radiation field was performed by 
R. B. Pietruszka [Ref. 1] using the dosimetry system which consists of 
Thermoluminescent Dosimeter (TLD) and associated TLD readers. The primary 
TLD's were СаҒо:Мп chip dosimeters and the TLD readers were Victoreen Model 
2600 and Model 2800M. This was the initial radiation mapping of the flash x—ray 


machine as a function of the system parameters. Since the Marx charge voltage is 


the primary variable controlled by the operator, this experiment was to measure 


exposure at various Marx charge voltages. 


When the radiation field was measured at the faceplate, one of the 
characteristics of the radiation pattern is that the maximum exposure did not occur 
at the geometric center of the anode assembly. The average radius of the relative 


maxima was 1.3 + 0.6 cm from the center. A typical pattern is shown in Figure 1. 
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Figure 1. Typical Pulsrad 112A Exposure Pattern at the Anode Face 
Plate. (Circles are at 100 kV Marx Charge, Stars at 75 kV Marx Charge) 


Variations in the exposures recorded were also observed between shots made at the 
same Marx charge voltage (see Appendix F). This variations are summarized in 
Table I for Marx charge voltage of 75 kV and 100 kV. Corresponding peak electron 
energy measured was 1.34 MeV for 75 kV and 1.66 MeV for 100 kV 


TABLEI: VARIATION OF EXPOSURE MEASURED AT GEOMETRIC 
CENTER OF ANODE ASSEMBLY [Ref. 1] 


Marx 

Voltage 15 KV 100 KV 
Peak Diode 

Voltage 1.34 MV 1.66 MV 
Minimum Peak 

Exposure 432 R 2040 R 
Maximum Peak 

Exposure 895 К 3090 К 
Меап 

Exposure 600 + 150 К 2600 + 340 R 


B. ITS/CYLTRAN 

ITS (Integrated TIGER series of coupled electron/photon Monte Carlo 
Transport Codes) is the most widely used computer code to predict radiation output 
generated from flash x—ray sources, as well as subsequent photon/electron transport 
and energy deposition in materials located downstream of the source. The release 
version 2.1 (February 1987) has been installed in the IBM main frame at Naval 
Postgraduate School and the release version 2.0 (July 1986) has also been installed 


on the micro VAX 3200 workstation in nuclear physics laboratory. 


The CYLTRAN code is a member code of the ITS system that is specially 
designed for applications having approximately azimuthal symmetry. CYLTRAN 
simulates the transport of particle trajectories through a three—dimensional 
multimaterial cylinder. For this thesis, only the CYLTRAN code was required. 
CYLTRAN, like all members of ITS system, contains the essential physics for 
generating and transporting the complete electron/photon cascade until these 
particles either are absorbed, escape from the system, or have their energies reduced 
below some user defined cutoffs. For more details of ITS/CYLTRAN see Appendix 
I 


C. PURPOSE 

In the present paper, the intensity of radiation field of the Model 112A 
Pulserad Pulsed X—Ray Generator is computed using the CYLTRAN code. The 
simulated outputs at the anode faceplate of the flash x—ray machine are compared 
with the experimental results in chapter III. The characteristic radiation pattern 
"dip" at the anode faceplate is utilized to investigate the electron source beam 
parameters such as effective cathode radius and angle of electron main flow. The 
radiation dosage in TLD's at the various axial distance is also calculated using 


CYLTRAN and compared with the design dose output. 


II. CALCULATION OF ABSORBED DOSE FROM MEASURED EXPOSURE 


А. PROCEDURES 

Since the CYLTRAN output is in terms of 'absorbed dose', the measured 
exposures should also be expressed as absorbed dose in order to compare each other. 
The absorbed dose in a medium can be determined from the measured exposures by 
following three equations (1) to (3). We have already the values of average exposure 
at 75 kV Marx charge and 100 kV Marx charge. In order to calculate the absorbed 
dose in a medium CaF:Mn, the following procedures were taken. 


(a) | Calculate the absorbed dose in air ( D...) using equation (1). 


s 
(b) Calculate the average photon energy absorption coefficients for air 
and CaF,:Mn. using equation (2) 


(c) Calculate the absorbed dose in CaF:Mn using equation I 


Under conditions of CPE. the absorbed dose is directly determined by the 


measurement of collision kerma, к. and exposure. X.[Ref. 5| 


(1) Des = ХӘ Eu D (for CPE) 
where B ы and D,.,. are in rads, and X in roentgens. 

The formal definition of CPE is "a condition that exists in a material under 
irradiation if the energies, number, and direction of charged particles induced by the 


radiation are constant throughout the volume of interest, Therefore, within such a 


volume, the sum of energies of all charged particles entering it is equal to the 
corresponding sum for all particles leaving it."[Ref. 12]. Moreover, if the same 
photon energy fluence V' is present in media A and B having two different average 
energy absorption coefficients (Hen /p) д апа (иеп /р) the ratio of the absorbed 


doses under CPE conditions in the two media will be given by 


(2) 2 о ы 
i % UP 


where 


( Hen/ P), and TAr ; average energy absorption coefficients 





Appendix D further elaborates on the terminology and definition used in these 


procedures. 


B. CALCULATION 

Since the values of exposure were measured the procedure (a) is simple to 
obtain. In order to calculate the average photon energy absorption coefficients, 
however the spectrum of photon energy fluence V'(E) is required. This distribution 
of photon energy can be obtained by running ITS According to ICRU(International 
Commission in the Radiation units and Measurements), the particle fluence of 
monoenergetic photons is defined as the differential quotient of the number of 


particles dN that cross a sphere of cross—sectional area da: 


¥ = dN/da 


This definition is not practical from a calculational point of view as it is based on 
the geometrical properties of a sphere. However, by multiplying both dN and da by 
the mean cord length l= 4V/A, a more practical definition can be obtained. As ldN 
is equal to the mean pathlength ds in the sphere, and dda is the volume dV of the 


sphere, the above equation becomes 


Y = ds / dV, 


which is independent of the shape of the volume, and can be interpreted as the 
"track length per unit volume". The distribution in photon energy fluence, V', can 
be calculated by ITS/CYLTRAN in terms of the "photon track length per unit 


volume" in a specific region of the test medium. 


~ 


There is another way of getting the incident photon energy fluence by ITS. 


The number of incident photons at a defined energy dE per unit area at (r,0), 


dN(E,0), can be calculated.(see Figure 2. Target Geometry). The incident photon 


fluence V' on the forward hemisphere of the photon target located at r is then 


Figure 2. 
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Let us take a spectrum of the photon energy transmitted at the anode face 
plate which is calculated by ITS/CYLTRAN. This photon energy spectrum is 
normalized to one incident electron and shown. in the Table IV and Figure 3. 
Obviously the average photon energy transmitted is less than the maximum energy 
and it is about 1/5 times of the maximum photon energy. More details of getting 
the overall photon energy spectrum will be discussed in the next chapter. The 
average energy absorption coefficients in air and in СаҒ2:Мп can be calculated from 
this photon energy fluence. This photon fluence is presented graphically in Figuree 
4. The absorbed dose in a test medium CaF:Mn from the measured values of 


exposure can also obtained using those average energy absorption coefficients. 


TABLE II:ENERGY SPECTRUM OF TRANSMITTED PHOTONS 
CALCULATED FROM CYLTRAN AT MARX CHARGE OF 100 KV AT 
ANODE FACE PLATE (normalized to one incident electron) 


Energy ( MeV ) Photons (#/SR—MeV) Photon 
Fluence(MeV/ Area) 

1.66 — 1.40 8.37 Е —05 1.26 Е 04 
140 — 1.3 9.37 E 04 12E 
ioe — 12 1.24 E —03 1.55 E 03 
12 — 1.1 3.90 E —03 4.03 Е —03 
mo — 1.0 3.68 E —03 3.86 E —03 
ШІ — 0.9 9:04 E, —03 9.96 E —03 
ШІ - 0.5 1.00 E —02 8.50 E —03 
Ше — 0.7 1.16 E —02 8.70 E —03 
0.7 — 0.6 1.81 E —02 ШЕ, 
ШІ - 0.5 3.17 E —02 1.74 Е —02 
ШӘ - 04 4.34 E —02 1.95 E —02 
0.4 — 0.3 1.15 E —02 2.50 Е —02 
US — 0.2 1.07 E —01 ОИЕ 502 
0.2 — 0.1 1.11 E —01 1.67 Е —02 
0.10 — 0.05 1.28 E —01 9.60 Е —03 
Total 4.5] E —01 ОЗЕ Ol 
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Anode Face Plate at a Marx Charge of 100 kV.(Normalized to Ore Incident 
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The result of calculation of the ratio of the two average energy absorption 


coefficients using equation (3) is as follows 
(%,77сарыма/ 0477. - 2:222 


If the measured exposure is X (roentgen), then the absorbed dose in СаЕ :Мп 15 
(0.876x2.222) x X (rad) for this particular photon energy fluence. The calculation of 
the average energy absorption is shown in the following Table V. This method of 
calculating the absorbed dose in the test material from the measured exposures will 
be applied in the following chapter and those absorbed doses will be compared with 


the simulation results. 


TABLE III: CALCULATION OF RATIO OF THE AVERAGE Ei eas 
ABSORPTIGN COEFFICIENTS 


Energy Photons Photon (en/p) ПЕТ ТТ! 
Fluence air aF4M n air CaF4Mn 
(MeV) (#/Sr—MeV) (MeV /area) (m^ /kg) (MeV /kg) 
1.661.4 8.31 E—5 1.26Е-4 0.0254 0.0247 3.20Е-6 3.11E-6 
Ко 9,3ТЕ-4 1.21Е-3 0.0266 0.0262 3.38E—5 . 39.90 B M 
| |2 1.24Е-3 1.55Е-3 0.0269 0.0265 4.17E-5 4Л11Е-5 
1.2 —1.1 3.90E—3 4.03E—3 0.0273 0.0267 1.10Е-4 1.05Е-4 
1.1 —1.0 3.68E—3 3.56Е-3 0.0278 0.0270 1.07Е-4 1.06Е-4 
1.0 —0.9 5.64E—3 ا‎ 0.0276 1.583E-4  1.48E-4 
0.9 —0.8 1.00E—2 5.50Е-3 0.0259 0.0251 2.46Е-4 2.39Е—4 
Е UT 1.16Е—2 о. (ОБ-ТІ 0.0255 2.94E-4 2.45Е-4 
0-7 —0.6 1.81 E—2 1.18E-2 0.0295 0.0290 3.45Е-4 . 3.13E-4 
0.6 —0.5 3.17E—2 1.74E—2 0.0295 0.0292 5.13E4 5.08Е—4 
у 0.4 4.34E—2 1.95Е—2 0.0296 0.0295 5.7тЕ—4 5.75Е—4 
0.4 —0.3 1.15Е-2 2,50Е-2 0.0290 0.0296 1.25Е-4 7.40Е—4 
учо 1.07Е—2 2.68E—3 0.0275 0.0306 1.37Е-5 8.20Е-5 
0.2-0.1 1.11Е—1 1.67E—2 0.0250 0.0397 4ЛТЕ-4 6.63Е-4 
0.1—0.05 кее 9.60E—3 0.0406 0.436 3.90E—4 — 4.19E—3 


TOTAL 1.36E—1 3.60E—3 ЭШЕ 


III. RADIATION DOSE ANALYSIS USING ITS/CYLTRAN 


The use of the Model 112A pulserad pulsed x—ray generator in radiation 
experiments requires knowledge of the characteristics of the output radiation. The 
experiment to measure the exposure at various axial/angular distance was done by 
Renee B. Pietruszka. Although the measured values have some variations from shot 
to shot, those measured exposures are the basis of this radiation dosage analysis. In 
order to compare the measured values with the results obtained by running the ITS, 
the experimental configurations at 75 kV and 100 kV Marx charge voltage were 
modeled with ITS/CYLTRAN codes. The dosimetry analysis using ITS/CYLTRAN 


was divided into three principal parts : 


i, Determination of the photon energy spectrum 


Ww 


Investigation of electron source beam parameters by comparison 
between measured radiation pattern and simulation results 
9. Determination of the intensity of the radiation at different axial 


distance. 


А. SIMULATION ENTRIES 
г Electron Accelerator Tube 
The electron accelerator tube is a field emission diode mounted in a 
vacuum. The cathode is a cylindrical stainless steel rod which tapers from a shank 
diameter of 3.2 cm (5/4 inch) to a tip diameter of 1.9 cm (3/4 inch). Either end of 


the cathode may be used for electron beam generation. The experiments to 
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measure exposures were made with the 3.2 cm (5/4 inch) end as the electron 
emission source. 

The anode, which acts as the bremsstrahlung target, is made from a sheet of 
0.381 mm (0.015 inch) tantalum. The tantalum is protected by a 0.64 cm (1/4 inch) 
thick aluminum faceplate. The anode-cathode spacing should have been 2.54 cm 
which is specified in the Pluserad 112A Operations Manual. However, when 
measured, the anode—cathode spacing was found to be 17.5 mm. The diagram of the 
electron accelerator tube is depicted in Appendix C. 

2 Electron Beam Energy 

a. Measured Diode Voltage and Current 
The Marx Output Voltage was monitored during the previous 

experiments using the two voltage monitors. These monitors (PIM —197) are CuSO, 
resistors The current across the diode was also monitored with the Model 199 
Fluxmeter. It consists of a loop that is mounted in the diode vacuum envelope. The 
changing magnetic field in the diode vacuum envelope induces a current in the loop. 
The resulting voltage pulse in the "B—probe" is output from the monitor. The 
magnetic field and the current in the diode can then be calculated. Table VI 


summarizes the results of data taken at Marx charging voltages of 75kV and 100kV. 


TABLE IV : MEASURED DIODE VOLTAGE AND CURRENT FOR 
PULSERAD 112A.[Ref. 1] 


Marx Charge Diode Voltage Diode Current 
(kV) (MV) (kA) 

Го 1:34 -= 0105 UO +1 

100 1:60. + 00: 2055 MI 
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The following two figures ( Fig. 5, Fig. 6 ) are photographs of typical diode 
voltage waveforms monitored by the Tektronix 7104 Oscilloscope. The half—width 
at half-maximum of both pulses is approximately.30 ns. These figures indicate that 
the electron beam is not monoenergetic. The energy distribution of the electron 
beam can be figured from these photographs. The photographs of diode current 
waveforms also were taken during the experiment. They are shown in Figure 7 and 
Figure 8. The half—width at half-maximum of these current pulses are 20 ns. The 
total electron flow during the pulse, which is needed for the calculation of absorbed 
dose in a test medium, can be extracted from these current waveform photographs. 
Calibration of these monitored pulses including monitor sensitivity and attenuation 


ratio of voltages in and out was described in Reference 1. 
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Figure 5. Photograph of Typical Diode Voltage Waveform at Marx 
Charge 75 kV (horizontal scale 20 ns/cm, peak voltage 1.2 MV). 
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Figure 6. "оюма of Typical Diode Voltage Waveform at Marx 
Charge 100 kV (horizontal scale 50 ns/cm, peak voltage 1.7 MV) 





Figure 7. Photograph of Typical Diode Current Waveform at the Marx 
Charge Voltage 75 kV (horizontal scale 20 ns/cm, peak current 11 kA)) 
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Figure 8. Photograph of Typical Diode Current Waveform at The Marx 
Charge Voltage 100 kV (horizontal scale 20 ns/cm, peak current 20.9 kA). 


b. Theoretical Diode Current and Voltage Outputs 
The voltage across the accelerator tube can be calculated by the 


following way.[Ref. 3| 


ЕЙ y. 
Mix ыу 


V» — 12x 1.2 x 0.9 x Vo 
where 
V, = Accelerator tube Voltage 
Va = Blumlein Pulse Charge Voltage 
Vo = Marx Charge Voltage 
Z; = Accelerator tube Impedance 
Zb = Blumlein Impedance 


The resulting diode current (1) can then be computed from Ohm's law: 


= ا 
= 


The following table summarizes the theoretical results of calculations for various 


Marx charging Voltages. 


TABLE V: THEORETICAL VOLTAGE AND CURRENT PARAMETERS 
FOR PULSERAD 112A [Ref. 1] 


Vo (Marx) Vp(Blumlein) Vi (Tube) It Сезне 
(kV (MV) (MV) (КА) 

(Б; 0.97 1.04 20.90 
100 1.30 1.39 2s 


There is a significant difference between the calculated and measured values 
(see Table VIII). The measured impedance (Zt = Vt/It) was found to be 122 ohms 
at Marx Charge 75 kV and 80 ohms at Marx Charge 100 kV which are different 
from the specified electron tube impedance value of 50 ohms. This may be caused 
fron the fact that the cathode—anode seperation was smaller. The further 
investigation of diode spacing is needed. The measured values were used for the 
input data to analyze the dosage patterns in the following Chapters. The next table 


campares the theoretical voltage/current with the measured voltage/current. 


TABLE VI: COMPARISON BETWEEN THEORETICAL DIODE 
VOLTAGE/CURRENT AND MEASURED DIODE VOLTAGE/CURRENT 


Marx Voltage Diode Voltage Diode Current 
(kV) Measured Calculated Measured Calculated 
(MV) (MeV) (kA) (kA) 
15 1.34+.03 1.04 I 0l 20.9 
100 1.66+.03 1.39 20 5t] 27.9 
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B. DETERMINATION OF PHOTON ENERGY SPECTRUM 

The first step to obtain the radiation output of the Model 112A Pulserad 
Pulsed X—Ray Generator is to determine the Photon energy spectrum. The Photon 
energy spectrum was obtained by running ITS. The measured parameters of the 
electron beam accelerator and its measured beam energy, which are specified in the 
previous section A, were used as input data. Although the input parameters should 
be the same as in the Flash X—ray machine, there are several unknown parameters 
that are necessary for the simulation input. 

The first unknown parameter is the radius of the electron source beam. The 
beam radius may not be the same as the radius of the cathode. In other words the 
effective cathode radius may be larger than the physical cathode radius. Futher 
investigation on this point will be done in the next section B. Three sets of the 
beam radii were entered for the overall photon energy spectrum which is normalized 
to one incident electron. 

The second parameter is the beam propagation angle. Because the dominant 
electron flow mode may be parallel to the electrostatic field lines of force, the source 
beam would be truncated at a certain angle with respect to the center direction. The 
key words "isotropic" and "cosine-law" are used in ITS/CYLTRAN for the source 
electron angular distribution. The word "isotropic" defines angular distribution of 
source particles as being isotropic with respect to the reference direction, and the 
word "cosine-law" defines angular distribution of source electrons as being 
proportional to the cosine of the specified angle with respect to the reference 
direction. The distribution is truncated at a given angle in both kevwords. The 


angles 30, 45 and 60 degrees for both "isotropic" and "cosine-law" were entered to 
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investigate the overall photon energy spectrum. Results for monodirectional flow 
with all electrons parallel from cathode to anode is also calculated. 

The third ambiguous factor is the energy distribution of the electron beam. 
Photographs of the diode voltage wave form were taken for the duration of the 
pulse. The measured diode voltage and current are different from the theoretical 
values. Because the photographs of diode current waveform and diode voltage 
waveform are available, the measured values of diode voltage were used for the 
overall photon energy spectrum. From those photographs the energy spectrum of 
electron source beam is simplified for the simulation input by applying a histogram 


fit to the typical diode waveform( see Fig. 9). 
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Figure 9. Distribution of the Beam Energy Content for 75 kV 
Marx Charge(upper) and 100 kV Marx Charge(lower). 
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l. Input 
The radius of the electron source beam, the angle of the electron flow 
mode, and the beam energy distribution used for determining the overall spectrum 
of transmitted photon energy were shown in the Table VIIII. The simulation 
geometry is depicted in Figure 10. The actual input file for running ITS/CYLTRAN 
is in Appendix A. 


TABLE VII : INPUT DATA 


beam radius — 5/8inch(1.86 cm) 
— ] inch(2.54 cm) 
— 3/2inch(3.81 ст) 


source electrons 
diode space ШІ сіп 
anode TA : 0.0381 cm thick 


AL : 0.635 cm thick 
anode radius 10 cm 


beam energv 
— monoenergetic beam 
1.34 MeV at Marx charge 75 kV 
1.66 MeV at Marx charge 100 kV 
— Spectrum type beam 


see Fig.8 
cutoffs 0.05 MeV 
angular 
distribution — monodirectional(perpendicular to anode) 
— isotropic 30, 45 and 60 degrees 
— cosine—law 30, 45 and 60 degrees 
histories 35,000 





Figure 10. Geometry Configuration entered into CYLTRAN 
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2. Output 

The outputs of the transmitted photon energy spectra are normalized 
to one incident particle. The photon energy spectrum for both the monoenergetic 
electron source beam and the spectrum type source beam are presented in the Fig 
11. and Fig 12. The photon energy spectrum of the spectrum source beam is about 
twice that of the monoenergetic beam at both 75 kV Marx voltage and 100 kV marx 
charge. This result is reasonable as the average energy of the spectrum source is 
much less than that of the monoenergetic beam. 

Although the electron source beam is expanding in the vacuum tube with a 
certain angle with respect to the reference direction, the overall energy spectrum of 
the transmitted photon remains almost the same. The comparison of the photon 
energy spectrum for the different angular distribution of electron beam is shown in 
figures 13, 14 and 15. These figures shows that the overall transmitted photon 
energy spectrum is independent on the angular distribution of electron source beam. 

For three different cathode radii the photon energy spectra also appear to be 


the same. This phenomenon is presented in Figure 16. 
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electron Beam and Spectrum Source at a Marx Charge of 100 kV 
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Figure 13. Transmitted Photon Energy Spectra for Spectrum Source 


Beam Having Different Angular Distribution(cosine law) at a Mmarx Charge of 75 
kV 


28 


BE 





3 
<Я 
> 
m 
D 
ші 

S 3 5 5 5 

о e о O O 

S{HUOAM/AINXKYS7IIQUNN 
Figure 14. Transmitted Photon Energy Spectra for Spectrum Source 


Having Different Angular Distribution(isotropic) at a Marx Charge of 75 kV 


29 


E 





2 
© 
>” 
c 
Ф 
c 
Ш 

a 

о 

о 

e[onJDg/ASym*HS/19quanw 
Figure 15. Transmitted Photon Energy Spectra for Spectrum Source 


Beam Having Cosine-law distribution and Isotropic Disrrbution at the Marx 
Charge of 75 kV. 
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Figure 16. Transmitted Photon Energy Spectra for Different Cathode 


Radii at 75 kV Marx Charge (spect1; 1.8575 cm, spect2; 2.54cm, spect3; 3.81 сіп) 
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C. INVESTIGATION OF THE ELECTRON SOURCE BEAM 
PARAMETERS 
For the further analysis of the radiation output of the Model 112A x-ray 
generator, the unknown parameters of the electron source beam should be specified. 
Since the ITS/CYLTRAN codes can calculate the absorbed dose in the test material 
CaF:Mn, the radiation dosage patterns at the anode faceplate can be extracted. By 
comparing the measured exposure pattern with the simulated dosage patterns of 
ITS/CYLTRAN output, electron source beam parameters which are difficult to 
obtain by direct measurement can be determined. 
T CYLTRAN Input Parameters 
The following three beam parameters are the major factors to be 
specified. 
(a). the electron source beam radius(effective cathode radius) 
(b). the angle of the dominant electron flow 
(C). the peak energy of the electron source beam(or the peak 
current of the electron flow) 
All the other input data are the same as those entered in previous section A (see 
Table VIIII) and various values of the above three parameters are entered. The 
modeled configuration of the electron accelerator tube for the calculation of the 
radiation pattern in the test material is depicted in Figure 17. The TLD is divided 
into 10 spherical shell "zones" to have the radial dosage pattern. In order to get the 
closest simulation results a 'trial and error' method was applied in entering the 


input parameters. 


If the absorbed dose in air, D can be calculated directly by 


аш” 
CYLTRAN, then the comparison is easy and the calculated exposure 1s obtained by 
the simple relation (1) in chapter II. Within a' single Monte Carlo calculation, 
however, it is difficult to include both the generation of photons in the target and 
subsequent energy deposition from those photons in the small volume of air. The 
difficulty is caused by the low photon conversion efficiency for the electrons in the 
energy range considered here and by the low probability for photon interaction in 
small volume of air [Ref. 13]. Since this is also true for the test material СаҒо:Мп, 
the thickness of TLD was enlarged for the simulation input. The actual thickness of 


TLD in the experiment was 0.035 inch (0.089 cm). However, a thickness of 0.5 cm 


CaF5:Mn was modeled in CYLTRAN. 
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Figure 17. Geometry Configuration entered into CYLTRAN for 


Calculation of Absorbed Dose in TLD's. 
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D 75 kV Marx Charge Series 
а. Development of Dip at the Faceplate 

One of the unexplained characteristics of the measured 
exposure pattern is that it has a dip at the geometric center of the anode faceplate. 
The simulated output of the ITS/CYLTRAN also has the same pattern at a specific 
angle of the dominant electron flow. For a Marx charge of 75 kV, Figure 18 shows 
the absorbed dose in TLD normalized to one incident electron at different source 
beam flow angles using the cosine-law distribution. The radius of the cathode in 
this case was 2.54 cm. At cosine-law 30, 45 and 60 degrees the dip appears at center 
and the pattern is similar to the measured exposure pattern. Cosine—law 45 degree 
simulation shows the most extreme dip pattern. The source beam distribution of 
cosine—law 60 was the most probable configuration when compared to the 
normalized patterns. The calculation was normalized, as will be explained in the 
section B, by dividing the calculated and measured absorbed dose by the maximum 
value. For various beam radii of the cosine-law 60.0 degree distributions, the dosage 
patterns were compared with each other in Figure 19. The cases of isotropic 
distribution of the electron source beam were also simulated and the characteristic 
center 'dip' did not show up (see Figure 20). In Figures 18, 19 and 20, the 
integration of the dosage throughout the whole distance is not constant. The dosage 
patterns in the anode (Ta converter and Al absorber) do not show a dip, but the 
characteristic dip appears in the TLD's chip. Figure 21 presents the development of 


a dip through the anode. 
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Figure 18. Radiation Dosage Pattern in TLD's for Different Source Beam 


Distribution (cosine-law) with a Marx Charge of 75 kV(normalized to one incident 
electron) The beam radius is 2.54 сіп 
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Figure 20. < Radiation Dosage Pattern in TLD's for Different Source Beam 


Distribution (isotropic) with a Marx Charge of 75 kV(normalized to one incident 
electron) The beam radius is 2.54 cm 


38 


6.0E—- 002 


2:05 002 


dose(Grad /Coul) 


0.0E+0C0 


О 
distance(cm) 


SUE OS 


dose(Grad/c) 


0.0€+ 000 
9 
DISTANCE (CM) 





Figure 21. ۰ Radiation Dosage Pattern in TA(upper) and in AL and TLD 
(lower) for the Cosine-law 60 degree Beam Distribution with a Marx Charge 
of 75 kV (normalized to one incident electron) 
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b. Output and Comparison 

In order to compare the simulated absorbed dose which is 
directly obtained by running CYLTRAN with the measured exposure which is 
obtained from the experiment, the exposure values measured should be converted to 
absorbed dose values in the test material. The procedures for calculation of the 
absorbed dose from the measured exposure were followed as explained in Chapter II. 
The photon energy fluence in each zone obtained by CYLTRAN in terms of "photon 
flux distribution (track—length/Vol—MeV—particle)" at 75 kV Marx charge is 
presented in Table XI. Since the absorbed dose in TLD calculated by CYLTRAN is 
normalized to one incident electron, these computed absorbed dose also have to be 
converted to absorbed dose for the each electron propagating during the pulse. From 
the current wave form taken during the experiment (see Figure 7 and Figure 8), the 
number of electrons can be obtained by integrating the current for the duration of 
pulse. 

Large difference between the absorbed dose converted from the 
experimental exposures averaged and the absorbed dose calculated by CYLTRAN 
were observed. This comparison 1s shown in Figure 22. The calculated absorbed dose 
is factor of 2 to 5 times smaller than the measured absorbed dose. There may be 
some experimental reasons for this large difference. One possible reason would be 
that the operation time was more than the 20 nsec specified in the operations and 
maintenance manual. Since the system 1s based upon time-dependent electrical 
breakdown, the longer switch time may result in higher radiation output. Another 
error source may be the unshielded TLD's used during the experiments. Unshielded 


dosimeters do not provide electron equilibrium in a high energy gamma field. The 
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ratio of the indicated exposure and actual exposure in pressed chip TLD's ranges 
from 2 to 10 [Ref. 7]. So unshielded TLD's may result in higher exposures than 
actual exposures. For this reason the measurements of Reference 1 should be 
repeated using shielded doimeters. 

Because there was such significant difference between the 
CYLTRAN results and the measured absorbed dose, a normalized dosage pattern 
was produced. This normalization was performed by dividing each absorbed dose by 
the maximum values so that the dosage patterns could be compared. The 
consequences of this normalization are depicted graphically in Figure 23, which 
compares the experimental results and CYLTRAN outputs. 

From these comparisons the possible source beam parameters, 
which were entered and produced the closest radiation dosage pattern to the 


experimental results, are summarized in the following table. 


TABLE VIII : SOURCE BEAM PARAMETERS 
(produced the most similar dosage pattern to experimental result at the Marx 
charge of 75 kV) 


Simulation Input 


Beam Radius 294 cm 
Electron Flow cosine-law 609 
source Energy peak energy 
spectrum 1.34 MeV 
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3. 100 kV Marx Charge Series 
The same method of analysis and calculation procedures were used for 
the Marx charge of 100 kV to investigate the electron source beam parameters. The 
simulation input parameters given in table XII produced the most similar dosage 
pattern to the experimenatal results. The comparison between the CYLTRAN 
output and experimental result is presented in Figure 24. There is also a huge 
difference between them. The same normalization to the 100 kV series was done and 


compared with each other and presented in Figure 25. 


TABLE X : SOURCE BEAM PARAMETERS 
(produced the most similar dosage pattern to experimental result at the Marx 
charge of 100 kV) 


Simulation Input 


bean radius оо 
electron flow cosine-law 459 
source energy peak energy 
spectrum 1.66 MeV 
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Figure 24. ~ Comparison of the Absorbed Dose Between the CYLTRAN 
Outputs and Experimental Result at a Marx Charge of 100 kV 
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D. PREDICTION OF THE RADIATION DOSAGE 
The design specifications of the Model 112A Pulserad Pulsed X-Ray 
Generator is provided in the following table as listed in the Operations and 


Maintenance Manual. 


TABLE XI: SYSTEM SPECIBICATIONS 


Marx Charge Voltage 100 kV 
Output Voltage 1.7 МУ 
Pulse Width 20 ns 
Exposure at 0.5 m from anode 8 roentgen 


The same simulation input parameters of the electron source beam, as given 
in Table XII, were entered. The simulation model geometry is depicted in Figure 26. 
The AL equilibrator of 0.2 cm thickness was placed to establish Charged Particle 
Equilibrium (CPE). This equilibrium thickness is proper for our photon energy 
transmitted [Ref. 9]. Figure 27 shows the radiation pattern in TLD's in rad at the 
axial distance 0.5 m fron the anode assembly. The procedures in Chapter III were 
followed to convert the design dose specification in roentgen to the absorbed dose in 
rad in TLD's. in order to compare this design dose output with the CYLTRAN 
result. The comparison between the CYLTRAN output and the design dose output 


is presented in Table XIV, showing that they are in good agreement each other. 
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TABLE XII: DOSE COMPARISON BETWEEN CYLTRAN OUTPUT AND 
DESIGN SPECIFICATION (at the axial distance 50 cm from anode) 


CYLTRAN output Design dose 
8.5 rad 8.2 rad 
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Figure 26. ` Geometry Configuration of CYLTRAN for Calculation of 
Absorbed Dose in TLD's at the axial distance 50 cm from anode 
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IV. CONCLUSIONS AND RECOMENDATIONS 


The CYLTRAN simulation modeled for the Model 112A Pulsrad Pulsed 


X-ray Generator can be concluded as follows: 


p The overall photon energy spectrum transmitted from the x—ray 
machine is independent of the electron source beam parameters such 
as beam radius and electron angular distribution. 

2 The source beam was found to have angular distribution of electrons, 
confined within 60 degrees and proportonal to the cosine of the angle 
with respect to the cathode center axis at a Marx charge of 75 kV 
and 45 degrees at 100 kV Marx charge. 

3k The CYLTRAN calculations of the absorbed dose in TLD's were 
considerably smaller than the measured values at the anode faceplate. 

4. However, the CYLTRAN output agreed with the design dose 
specification at the distance 50 cm from anode. This agreement may 


be an indication that the same code was used by the manufacturer. 


To make effective use of the Pulserad 112A a experiment of measurement of 
the radiation dose with establishment of Charged Particle Equilibrium at the 
various axial distance is recommended. The experimental method is described in 


detail in Reference 9. 
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APPENDIX А. SAMPLE INPUT FILES 


sample input file to generate a cross section tape. 


Energy 1.7 

Material Ta 

Material Al 

Material Ca 0.495 F 0.484 Mn 0.021 Density 3.18 
title 


1.7 MeV Cross Section for Ta, Al and CaF,:Mn 
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Sample Input File to excute ITS 


Echo 1 
Title 
Mat; Ta, Al, СаҒо:Мп His; 40000 Energy; 1.34(spectrum) 


Spectrum 7 


1.0 0.643 0.439 0.246 0.096 0.032 0.0 
1.34 1.0 080.604 0 21070 
Cutoffs 0.05 0.05 
Position 0.0 0.0 0.0 
Radius 2.54 
Direction 0.0 0.0 
Cosine-law 60.0 
Geometry 4 
0.0 1.75 00100 D0 
109 1.7881 0.0 100 ие 
1.7881 2423 00 100 2 TN 
2423] 279231] 2021 1 0 
Electron—escape 


Nbine 15 user 

1.41.3 1.2 1.11.0 0:9 0.8 0-7 0605 0.4 00 02 Oa 

Nbint 15 user 

1.2.3.5.1,10, 15.16. 20 30 25700 о ЧО И 
Photon—escape 

Nbine 15 user 

1.4 1.3 1.2 1.1 1.0 0.9 0.80.71 000.5 01 03 0 

Nbint 15 user 

1. 2. 3. 5. 7. 10713. 1672039095 ODER тя 
Photon—flux 22? 31 

Nbine 15 user 

1.4 1.31.2 1.1 10 0.90.8 0706050070 OTEC 
Pulse—height 22 

М рше 

Tool УЕ 
Histories 40000 


31 
user 
21.1 1.0 0.9 0.8 0.7 0.6 Ula G24 Ora RN 


94 


APPENDIX B 


MASS ENERGY-ABSORPTION COEFFICIENTS р, /p (cm^/g) FOR 
VARIOUS MEDIA [Ref. 5] 


Photon Energy Air CaF,:Mn 
0.01 4.61 ЭТЕТ 
0.015 1.27 16.1 
0.02 0.511 6.86 
0.03 0.148 2.03 
0.04 0.0668 0.850 
0.05 0.0406 0.436 
0.06 0.0305 0.256 
0.08 0.0243 0.118 
0.10 0.0234 0.0697 
0.15 0.0250 0.0397 
0.2 0.0268 0.0317 
0.3 0.0287 0.0296 
0.4 0.0295 0.0295 
0.5 0.0296 0.0293 
0.6 0.0295 0.0290 
0.8 0.0289 0.0281 
1.0 MeV 0.0278 0.0270 
1.5 MeV 0.0254 0.0247 
2.0 MeV 0.0234 0.0229 


Of importance in this table is the significant dependence on the low energy 


component. The cutoff energy of photon should be as low as possible in running the 


ITS/CYLTRAN codes. 
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APPENDIX C 


DIAGRAM OF THE PULSERAD 112A ELECTRON 
ACCELERATOR TUBE 





а= 1 50 ср 
b = ІШІСІ Сі 
с = 0005 Cm 
d = 0635 CT 
eol ccm 
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APPENDIX D. TERMINOLOGY AND DEFINITION 


А. EXPOSURE AND KERMA 
Exposure is symbolized by X, and is defined as " the ionization equivalent of 
the collision kerma К, in air, for x-ray and 4-rays." For monoenergetic photons 


К, is related to the energy fluence Y by the mass energy—absorption coefficient (i A 


/ P)g 7 


и 
(1) ГЕ ТЕ: 





where 
V — energy fluence (J/ cm?) 
Hon = linear energy absorption coefficient (1/cm) 
p = density (g/cm) 


We can write that the exposure at a point due to an energy fluence ¥ of 


monoenergetic photons of energy E is given by 


(2) x =v (E 


jm ER кез Ш-СЕС,,/ 350 
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where 
X is the exposure in C/kg 
V is expressed in J/ m* 
2 2 
(enl P) E air is in m^/kg 
(e/ МЕ ИИ ОКЕ 


The roentgen (R) is customary unit of exposure. It is converted by 


—4 


X (C/kg) 2 2.580 « 10 ^ X (R) 


X (R) = 3876 X (C/kg) 


If a spectrum of photon energies W'(E) (in j/m?—keV ) is present at the point of 
interest P, and if (527 ЭБ ET (in m^ /kg) is the energy—absorption coefficient as a 


function of photon energy E for air, then the exposure at P will be given bv 


e/ W),.. V'(E) dE 


Emax 
(3) x | 


Е-0 (Hon! PB air ( air 
where 

(cee = ООА 

V'(E) — distribution in photon energy fluence; J/(cn^—keV) 


dE is in keV 
X is in C/kg 


В. ABSORBED DOSE AND KERMA 
The absorbed dose (D) is the quotient of de by dm, where de is the mean 
energy imparted by ionizing radiation to matter of mass dm. The unit for absorbed 
dose is the gray( Gy). 
кила укр =н су 
The absorbed dose represents the energy per unit mass that remains in the material 
The establishment of " charged particle equilibrium (CPE) " enables us to 
relate " absorbed dose " to " kerma". Under conditions of CPE, the absorbed dose is 


directly determined by the measurement of collision kerma. Ref. 5| 
(4) D, —(K),,-0876X (forCPE) 


where (K and D,,.. are in rads, and X in roentgens. 


pu 

Moreover, if the same photon energy fluence V is present in media A and B 
having two different average energy absorption coefficients (jen /p) қ and (Hen lP) 
the ratio of the absorbed doses under CPE conditions in the two media will be given 


by 


(5) EE gi — = (for CPE) 


where 


(Hen/P) and (ҮГИ ; average energy absorption coefficients 


C. CHARGED PARTICLE EQUILIBRIUM 

The formal definition of CPE is "a condition that exists in a material under 
irradiation if the energies, number, and direction of charged particles induced by the 
radiation are constant throughout the volume of interest, Therefore, within such a 
volume, the sum of energies of all charged particles entering it is equal to the 
corresponding sum for all particles leaving it." 

To understand the charged particle equilibrium (CPE), consider a small 
volume of material, a mass of Am. Let this volume of material be traversed by a 


beam of ionizing radiation. The absorbed dose in the material is by definition 


AE 


AD — D 
Am 
where 
AED = AE - AE, ; the energy deposited in the material 
АЕр, АЕр ; the energy entering and leaving the material during 


some time interval 


In the case of incident photon radiation only the energy transfer to matter is a two 
step process: 

1. Photons impart energy to electrons via the photoelectric. compton 
scattering, and pair production processes. 

2. Electrons impart energy to matter via excitation, ionization, and 


elastic scattering processes, 
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Since the energy deposition process involves only photons and electrons, the energy 


deposited can be redefined as 
AE, = [AEp(7) + АЕр(е)] — [АЕ ( 7) *- ^E, (e)] 


where y and e refers to photons and electrons respectively. 
However, by definition, CPE exists when the energy carried into the mass by 


electrons equals the energy carried out of the mass by electrons. Hence, 
АЕр - АЕв(7) -АЕ4(7) (for CPE) 


The so-called equilibrium dose is then the ionization energy deposited in the 


material per unit mass 


Б АЕБ(1) = AE, (7) 


eq Am 


This expression for the equilibrium dose is equivalent to the definition of collision 


kerma.[Ref. 7] 
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APPENDIX E. ITS/CYLTRAN CODE SYSTEM 


The TIGER series of time-independent coupled electron—photon Monte 
Carlo transport codes is a group of multimaterial and multidimensional codes 
designed to provide a state-of—the—art description of the production and transport 
of the electron—photon cascade. All the codes are based primarily on the original 
ETRAN model which combines microscopic photon transport with a macroscopic 
random walk for electron transport. The basis of the ITS codes, ETRAN has been 


developed at the National Bureau of Standards. 


A. ETRAN MONTE CARLO METHODS 

The ETRAN model pertains only to schematization of the electron random 
walk that is used in the calculation. The methods used to generate electron 
trajectories go back to a paper published in 1963 [Ref. 8], and involve the sampling 
from relevant multiple scattering distributions. The code follows all generations of 
electrons and photons with energies up to 1 GeV and down to 1 keV in any target 
material. The result is a calculation which takes onto account primary electrons, 
positrons or photons,and all secondary radiations, including knock—on electrons 
from  electron-impact ionization events, electron bremsstrahlung. Compton 
electrons, photoelectrons, electron—positron pairs, annihilation radiation, and 
K-shell characteristic x-rays and Auger electrons resulting from electron and 


photon ionization events. 


1. Photon Transport 

Successive photon interactions are sampled individually in direct 
analogy to the physical process. Compton scattering, photoelectric absorption and 
pair production are included. The cross sections for these interactions are taken 
from the work of Hubbel done and have been organized into a database with which 
the necessary computation can be rapidly performed for any material. The photon 
Monte Carlo procedure is as follows. The distance to the next photon interaction is 
sampled from an exponential distribution governed by the attenuation coefficient. 
Then the type of interaction is then sampled from the appropriate relative 
probabilities. The history of each photon is continued from collision to collision until 
the photon either absorbed, escapes the target, or its energy falls below a chosen 
cutoff value.[Ref.6] 

2; Electron Transport 

The step pathlengths in ETRAN are chosen on two levels. First, 
major steps are defined such that, on the average, the kinetic energy of the electron 
is reduced by a constant factor. Then the major steps are divided into several 
equal—length sub—steps in order to keep the mean deflection angle small, which is 
required for applying the multiple—scattering theory. At the end of each sub-step 
the angular deflections are sampled. The net angular deflection from the combined 
effect of the elastic and inelastic collisions in a single sub—step is sampled from the 
Goudsmit—Saunderson multiple—scattering distribution. 

The energy loss due to multiple ionization and excitation collisions in 
a pathlength is described by the Landau distribution which is corrected by Blunk 
and Leisagang. In ETRAN, the collision energy loss for each major step is sampled 


from the Landau/Blunck-Leisegang distribution. The radiative energy loss for the 
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electron is determined by sampling the production of bremsstrahlung photons using 
a dataset of bremsstrahlung production cross sections developed by M. J. Berger 
and S. M. Seltzer. 

The production of knock-on electrons, whose energies are above a 
chosen cutoff value, is sampled according to the Moller cross section for collisions 
between free electrons. K—shell ionization events in each sub-step are sampled on 
the basis of Kolbenstvedt's cross section and either the emission of characteristic 
x-rays or Auger electrons are selected. For positron, most of the necessary 
cross—section information is available within the code, but separate sampling has 
not been implemented. 

Electron histories are terminated under a variety of conditions. A 
maximum cutoff value is selected so that no electron histories is terminated at a 
higher energy unless it escapes the target. Another energy is selected to be the 
minimum cutoff value to have no electron flowed to a lower energy. At energies 
between these two cutoffs, the history is terminated if the residual range is smaller 


than the distance to nearest boundary.[Ref.6] 


B. ITS CODE SYSTEM 

ITS is a powerful and user-friendly software package permitting accurate 
Monte Carlo Solution of linear time-integrated coupled electron/photon radiation 
transport problems,with or without the presence of macroscopic electric and 
magnetic fields of arbitrary spatial dependence. The ITS code svstem consists of 


four essential elements: 
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XDATA : The electron/photon cross section data file 


XGEN : The cross section generation program 
TS : The Monte Carlo Program file 
UPEML : A machine portable update emulator 


The heart of ITS is the Monte Carlo program file, which contains the eight 
codes of Table Il. Each of eight member codes will run on any of four 
machines-CRAY, CDC, VAX, or IBM. The update emulator program UPEML 
creates the various Monte Carlo codes for given system with any corrections to 
those codes that may be desired. Program XGEN generates the problem specific 
cross section data tape using file XDATA for referenced inputs and a user defined 
input file. The Monte Carlo codes then read in the cross section tape and process 


the user defined problem. 


TABLE XIII: ITS Member Codes 


Standard Enhanced Ionization Macroscopic Fields 
Codes and Relaxation 

( P Codes ) ( M Codes ) 
TIGER TIGERP 
CYLTRAN CYLTRANP CYLTRANM 
ACCEPT ACCEPTP ACCEPTM 
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One of the eight ITS code is CYLTRAN, which simulates the transport of 
particle trajectories through a three—dimensional cylinder. Only the CYLTRAN 
code was used for this project. The following steps are required to execute a 


CYLTRAN run. 


1. Generate a cross section tape based on the different type of materials 
contained in the cylindrical geometry of a problem 

2. Create an input file which lists all the input parameters required to 
calculate desired outputs. 

3. Submit the input file and the generated cross section tape to the ITS 


Monte Carlo codes to execute a run. 


In addition to certain diagnostic information, the default output consists of: 
і. Energy and number escape fractions (leakage) for electrons, 
unscattered photons and scattered photons. 
2: Charge and Energy deposition profiles. 


a An explicit statement of energy conservation. 
Theses data are sufficient to confirm the general partioning and conservation 
of charge and energy. In addition to the default output, a number of optional 


Outputs may be selected through the use of the appropriate keywords. These are: 


ІІ Escape fractions that are differential in energv for both electrons and 


scattered photons. 
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2; Escape fractions that are differential in angle for both electrons and 
scattered photons. 

d Coupled energy and angular distributions of escaping electrons and 
scattered photons. 

4. Volume-averaged energy distributions of electron and photon scalar 
flux for selected regions of the problem geometry. 

5. Pseudo—pulse—height distributions for selected regions of the problem 
geometry.( The key word pulse—height causes the spectrum of 


absorbed energy to be calculated for zones) 


Except for an initial diagnostic table containing accounting information on the 
various particle types, every output quantity is followed by a one or two digit 
integer that is an estimate of one—sigma statistical uncertainty of that quantity 


expressed in percent. 


C. COMPUTATION TIME 

The computation time depends on the choices of run parameters of the 
problem considered. The sample input file specified in the Appendix A has the 
incident histories of 40,000 and 1.66 MeV electrons follows down to 0.05 MeV 
through the 0.6731 cm thick anode assembly. This run involved about 880,000 steps 
for the primaries, and the sampling of approximately 29,000 knock—ons and about 
37,000 secondary histories. The run time required on the IBM mainframe and on the 


micro VAX 3200 workstation is summarized in the following Table III. The 


67 


computer run time depends largely on the secondary electron cutoff energy. 
Changing this cutoff to half would result in a run time smaller by about a factor of 


10. 


TABLE XIV : COMPUTATION TIME 
(for the input file in Appendix A) 


Major parameters IBM main frame micro VAX 


40000 histories 

Electron energy 
(1.66 MeV to 0.05 MeV) 

3 materials 56 minutes 45 minutes 
(Ta. Al, TLD) 
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Marx Charge 75 kV Series 


1000.0 


APPENDIX F. MEASURED EXPOSURE VARIATION 
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